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Abstract 

A systematic study of low mass dilepton production from tt — p interactions in a hot medium is 
presented. Applying finite temperature perturbation theory the dilepton rate, respectively the 
virtual photon rate, is computed up to order g 2 p . For dilepton masses below the p the two-body 
reactions tctt — ► p7*, np — > 717*, and the decay process p — > 7r7r7* give significant contributions. 
Non-equilibrium contributions to the thermal rate are estimated, including the modification of 
the particle distribution function with non-zero pion chemical potential. The comparison of 
the dilepton rate with the recent data measured in nucleus-nucleus collisions at SPS energy by 
the CERES Collaboration is also performed. It is shown that the additional thermal dileptons 
from 7r — p interactions can partially account for the access of the soft dilepton yield seen 
experimentally. 



1 Introduction 



Dilepton production is one of the interesting tools to study collective effects in strongly inter- 
acting matter produced in ultrarelativistic heavy ion collisions [p], @, 0, |J. This is particularly 
evident from the recent CERN experimental data at SPS collision energy. A significant en- 
hancement of the dilepton spectrum measured in A-A collisions as compared with p-p and p-A 
collisions is reported. Two experiments, CERES/NA45 || and HELIOS/3 |J have measured 
dileptons in the low mass range in S-Au, Pb-Au and S-W collisions respectively. Both these ex- 
periments have reported the excess of dileptons with invariant masses between 0.25 and 1 GeV 
in A-A nuclear collisions, compared to expected || or measured |J spectra in p-A collisions. 

These remarkable results indicate the appearance of collective phenomena in heavy ion 
collisions, like thermalization of the medium after the collision. Since the temperatures are 
not yet to be expected as high, in particular the reactions involving pions and rho mesons are 
important and have to be considered. In the kinematical window of low mass dileptons pion 
annihilation with the p in the intermediate state is a basic source of thermal dileptons, also 
responsible for the low mass dilepton enhancement reported by the CERN experiments. Indeed, 
recent theoretical calculations show that close to the p peak the thermal production rate due to 
pion annihilation is compatible with the experimental data |7|, ||. However, neither the size of 
the access below the p peak nor the shape of the distribution is quantitatively explained by this 
contribution alone. In a thermal medium, however, the partial restoration of chiral symmetry 
in hot and dense hadronic matter may modify the properties of vector mesons especially their 
mass or decay width ||, 10 1: in Ref. W\ it is shown that thermal dilepton production due to the 



pion annihilation process together with the assumption of a decreasing rho meson mass leads 
to a quantitative explanation of the enhancement of the low-mass dileptons observed by the 
CERES and by the HELIOS/3 experiments. 

The above results together with previous theoretical studies Jy], [12], |13| suggest the im- 



portance of pion scattering in a medium as a source of soft dileptons. In a thermal medium, 
however, 7r + 7r~ — > e + e~ is not the only process which has to be considered. For example the 
contribution from the two-body reaction 7r + 7r~ — > p7* — > pe + e~ does not have a kinematical 
threshold at 2m n , and therefore certainly dominates the basic pion annihilation process for 
dilepton masses M ~ 2m n . This example indicates, that a more complete analysis of the low 
mass dilepton spectrum originating from a possible thermal medium, requires further consid- 
erations. 

Previous calculations of the thermal rate of emission of direct, i.e. real, photons with energies 
less than 1 GeV have shown that there are two contributing reactions involving a neutral p: 
the annihilation process tt + tt~ — > p 0/ ~f and Compton like scattering 7r ± p° — > 7r ± 7 [TJ|. Both 



contribute to dilepton production, where the real photon is replaced by a virtual one. The 
importance of two-body processes for low mass dileptons has also been stressed in the case of 



production from the quark-gluon plasma [15 |. 



In this paper we compute the contributions to the dilepton production rate from n — p 
interactions for invariant masses up to the p peak. In sect. 2, after shortly treating the 7r + 7r~ — > 
7* Born rate, we discuss the "real" 2 — > 2 and p — > TTTf-f* reactions, at the two- loop level of the 
virtual photon selfenergy. Also the "virtual" two- loop corrections to the Born rate are derived. 
Sect. 3 extends the estimates of the dilepton rate taking into account out of equilibrium effects 
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by introducing the pion chemical potential. Finally, in sect. 4 we compare the dilepton rate 
with experimental data || in the low mass region. 



2 Dilepton production and tt — p interactions 

The thermal dilepton production rate is considered in lowest order of the electromagnetic 
coupling. First, the medium is taken at rest. The rate per unit space-time volume, dR = 



dN/d x, is related to the absorptive part of the photon self-energy tensor U^ v |2|, [LI], [16], [17 
For lepton pairs of invariant mass M, energy q$ and momentum q the relation is 

" Mqo)l^(q ,q), (2.1) 



dM 2 d 3 q/q 24vr 4 M 2 

where n is the Bose distribution function at temperature T = 1/(3 (cf. Appendix A). The 
tensor is constrained by current conservation q^H^q) = 0. 

We adopt the real-time formulation of finite temperature field theory to evaluate Imfl^. 



Using the Keldysh variant (with time path-parameter a = |L8|, [L9|, p0|1 ) the following useful 
relation holds in terms of the photon-selfenergy matrix 1T^ 6 , 

ImlT^go,^ = ^l-^mUqo,q). (2.2) 



2.1 Born rate 



In order to clarify the notation we give the one- loop expression for Il^ 2 . With the momentum 
labels indicated in Fig. 1 this pion-loop contribution with the temperature-dependent propa- 
gators (|A.1| ), summarized in Appendix A, reads 

mj 2 (q ,q) = -He) 2 / ^(p + p'f iD 12 (p)iD 21 (p') , (2.3) 

where p' = p — q. 




Figure 1: One-loop photon selfenergy diagram. The solid line denotes the pion 
Using the equilibrium distributions one recovers the standard expression for the Born term 



7r + 7r — > 7* — > l + I . For masses M near the p it is appropriate [|11] to implement VDM and 
to insert the pion electromagnetic formfactor 



m, 



(M 2 — m* 2 ) 2 + T 2 p m 2 p '' 
3 



(2.4) 



into the production rate. For our numerical analysis the parameters are chosen as m p = 
0.775 GeV, m* = 0.761 GeV and T p = 0.118 GeV in agreement with the measured pion 



electromagnetic form factor [II]. We do not consider the possibility of medium dependent 
resonance parameters, as done in 0. Also contributions due to the Al resonance are 
neglected. 

For soft dileptons with invariant masses M < m p the mass distribution per unit space time 
volume is obtained within a very good approximation from the rate with the heavy photon 
taken at rest with respect to the medium (also kept at rest) || |l0j 

« ^A***^ (2.5) 



dM 2 d A x y 1T' dM 2 d 3 q/q l 

where the Born rate for M ~ m p , with the Boltzmann approximation of the pion distribution, 
has the following simple form, 

| g -=o = ^|^(M)| 2 exp(-M/T). (2.6) 



dM 2 d 3 q/q ,q 96tH 

The more complete expression for the Born contribution with quantum statistics and arbitrary 
heavy photon momentum is, however, also well known in the literature (see e.g. [§]]). It will 
be used in the following numerical estimates. Throughout the paper we neglect the pion mass 
(except for the numerical analysis of the Born rate for M < m p ), and the masses of the leptons. 
Because of the many uncertainties we believe that these approximations are very well justified. 

Taking into account out of equilibrium effects in determining the rates, the expression ( |2.6| ) 
is multiplied by a factor 

1 + 25X, (2.7) 

with 5X = A — 1 in terms of the fugacity A (cf. Appendix A). This first approximation, i.e. 
staying near equilibrium for <5A <C 1, and neglecting more complicated dependences, should 
indicate the possible increase (A > 1) or the decrease (A < 1) of the dilepton rate due to 
non-equilibrium effects. 



2.2 Two- loop contributions 

The interaction of the charged pions with the neutral massive rho meson field p p and the 
electromagnetic potential A p is described by the Lagrangian [|14] 

L = \D^\ 2 - l Pia ,pr + \m 2 pPvP v - i V (2.8) 

where D p = d p — ieA p — ig p p p is the covariant derivative, $ is the complex pion field, p pv is the 
rho and F pv is the photon field strength tensor. As is well known, the pnir coupling is rather 
large g p /4:7T ~ 2.9; nevertheless we attempt an effective perturbative treatment up to two-loops. 
However, no resummation of high temperature effects, e.g. comparable to the hard thermal 
loop expansion for QCD |22| , ^3], is implemented in our approach, which maybe compared with 



an analogous fixed-order calculation for dilepton production from a quark-gluon plasma |24j . 



From the Lagrangian ( |2.8j ) and using the closed-time-path formalism we calculate the 

dilepton production rate produced in a thermal pionic medium at the two-loop level. Typical 
diagrams are shown in Figs. 2 and 3. The two tadpole diagrams of 0(g p ) are not included, 
since they do not contribute to the discontinuity of the photon selfenergy (|2.2|). 
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2.2.1 Real p° processes 

The processes involving real p°'s, namely tttt — > /ry*, 7rp — > 717* and p — > 7T7T7*, are expected 
to be important for dilepton masses M below the p-peak. These contributions are obtained by 
cutting the two-loop diagrams (Figs. 2 and 3) such that the p is put on-mass shell. 



P + k 




P' 

Figure 2: Two-loop diagram: self energy insertion. The labels a,b =1,2 denote the type of np 
vertex. The solid line with momentum label k corresponds to the p 





Figure 3: Two-loop diagrams: vertex type corrections 

We start to discuss explicitly the contribution from the diagram shown in Fig. 2, where the 
selfenergy correction 5D 12 (p) is inserted into the pionline with momentum p. The expression 
for 5D\ 2 {p) is derived in Appendix B using the Dyson equation. 

In this section we assume (thermal and chemical) equilibrium distributions ( |A.3| ), and con- 
sequently we use 5DH(p) given in (P-5| ). The discussion of out off equilibrium effects in the 7* 
rate is postponed to sect. 3. 

Using the second term of ( |B.5| ) we calculate the "real" correction to (|2.3|) 



^ S V,?1 = -He) 2 /-^(p + p0 2 p(^) (-iIl 12 (p))iD 



2i(p'), (2.9) 



where the pion selfenergy iU 12 (p) at one-loop order due to irp interactions (|2.8| ) is given by (cf. 
|BT6|and|C^), 

^ / U U \ 

(2.10) 



iU 12 (p)=g 2 p J^L{2p + kY (-9~ + ^rj (2p + kYD 12 (p + k)D 21 (k). 



We note that in Q2.9Q we only keep the real p° contributions. 

Inserting the propagators ( A.l ) with (|A.9| ) it is convenient to separate the different kine- 
matical configurations, which appear due to the presence of the sign function in ( |A.9| ). We 
consider in more detail the process of p-emission, i.e. vr(p + ) + 7r(p_) — > p(k) + J*(q), where 
we relabel to momenta of Fig. 2 as: p' — > — p_, p — > p+ — k, and introduce (positive) energies 



k° 



E P , p°_ 



Pa 



+ 4eV 



E + . We find the following contribution to i8Ii 



n(E, 



(2tt) 3 2£+ v +/ 7 (2tt) 3 2£_ 
(2tt) V(p + +p_-k- g)(2p_ -q) 2 P( 



f2vr) 3 2£„ 



P 




(2.11) 



x(P + -k)\ 

expressed in a form familiar from kinetic theory ||: it consists of the square of the matrix 
element for the two-body process mi — > fry*, here to one of the exchange diagrams plotted in 
Fig. 4, with the integration over the phase-space of the participating pions and the p-meson 
properly weighted by the thermal distribution functions n(E). Next we also integrate with 
respect to the photon momentum q as described in Appendix D. 



71 



Y 



* K 



K 




Figure 4: Matrix elements for the process rnr — > /ry* arising from cutting the two-loop diagrams 



An analogous, straightforward treatment gives the corresponding contributions to the real 
p° processes from the two-loop vertex type diagrams of Fig. 3. Finally we sum all these contri- 
butions, including those obtained by interchanging momenta, like the selfenergy contribution 
of ( |2.9| ), but interchanging p and p' in Fig. 2. 

As a check the same final result is obtained when the square of the sum of the three 
matrixelements illustrated in Fig. 4 is taken. This approach has e.g. the advantage that 
current conservation is rather easily checked on the level of these matrixelements. 

Here we quote the result. Mainly for simplicity we use the Boltzmann approximation 
which amounts to underestimate the dilepton rate: in this limit we replace e.g. in ( |2.11| ) 



n(E + )n(E^) — > exp(—E + — EJ)/T, and neglect stimulated emission, l + n(E p ) ~ 1. The dilep- 
ton rate due to the real p° processes is best expressed using the Mandelstam variables defined 
in ( |D.4|) ; together with the kinematical approximations discussed in Appendix D, we obtain for 
the dilepton mass range M < m p , 



dM 2 d i x ~ 24vr 4 M 2 \ 



ttT 3 
2m 3 



x 



dse 2m " T + / dse~ mp/T 

(m p +M) 2 JO 

m 2 M 2 /l 1\ (m 2 n + M 2 ) 2 + m 2 M 2 /2 



du + 2 P / dte 2m P T 



tu 



{ml + M 2 ) 



1 1 

t u 
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where 



« + =\{m 2 p + M 2 - s) (-) i ^/( a -(m p + M)2)( a -(m p -M)2), 
u max = m 2 M 2 /t, u min = m 2 + M 2 - t, (2.13) 
with s + t + u = m 2 + M 2 . The ^-integration may still be performed analytically. 



The kinematic boundaries of integration given for the processes, which are summed in ( 2.12T) 



in the order of tttt — > p7*, p — > mvy*, irp — > n-y*, are illustrated in the Mandelstam plot of Fig. 8. 
The result contained in fl2.12j ), for the squared matrixelements of Fig. 4, is sucessfully 



compared with the one derived in [[TJj] for the case of the two-body processes involving a real 
photon. 

It is important to note that due to the non- vanishing p mass and the heavy photon, M > 0, 
no mass singularities appear in ( |2.12| ), even when setting the T = pion mass to zero. However, 
for the evaluation of the 7rp — > 7T7* contribution, the principial value prescription leading to 
a well defined integral, e.g. in ( |2.12j ) P(l/u) = lim £ ^ u/{u 2 + e 2 ), and correspondingly for 
P(l/w 2 ) [0, has to be taken into account, in order to correctly treat the behaviour near u ~ 0, 
which is covered by the kinematical domain of this process (Fig. 8). For the processes tttt — > p7* 
and p — > 7T7T7* the P prescription is not explicitly required. 



We plot the rate (|2.12| ) below the p-peak in Fig. 5 as a function of the dilepton mass M for 
M > 0.2 GeV at fixed temperature T = 150 MeV. For comparison the Born rate tttt — > 7* is also 
included. We find that, as expected, the two-body reactions tttt — > p7* and irp — > ttj*, together 
with the decay channel p — > tttt^*, dominate the dilepton rate for low masses^. Thus, when 
discussing the thermal yield in the context of recent experimental data on dilepton production 
in heavy ion collisions, it is necessary to include the above higher order processes. 

It is worth to mention that admitting in-medium effects on the p meson mass or width, 
as recently proposed in Refs. [7, 9, 10], could substantially influence the overall thermal rate 
(solid curve) shown in Fig. 5. In particular, assuming a decrease of the p meson mass by only 
100 MeV would imply an increase of the thermal rate (solid curve of Fig. 5) by a factor of 2. 



2.2.2 Virtual p° processes 

For heavy photon production at 0(g 2 ) we have, in order to correct the Born rate ( |2.6| ), to 
include virtual contributions, which arise from the processes shown in Figs. 2 and 3 by cutting 
the diagrams in the proper way only through pion lines, without cutting the p line. For the 
production of real photons |L4] these contributions are absent. 



In some detail we describe our estimate, including the approximations, for the selfenergy 
diagram (Fig. 2). Here we take the first term of ( [B.5| ) and evaluate (cf. 

*5H <T rt5 V, q) = He) 2 J0^(P + P'f e(Po) n(p ) 5'(p 2 ) ReH(p) iD 21 (p'), (2.14) 



^^We note that the curves do not strongly change when the kinematical approximations performed in Appendix 
D are relaxed, e.g. for 7T7t — > p^j* at most a 10% change is numerically found when the p is not kincmatically 
kept fixed at rest 
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Figure 5: Dilepton rate as a function of M at fixed T = 150 MeV. The solid curve represents 
the sum of the real and the virtual processes 

where now the real part of the pion selfenergy H(p) enters, and equilibrium distributions are 
used. 

The temperature-independent part of Rell is absorbed in the definition of the T = pion 
mass (which we take approximately as m n = 0). The temperature dependent part in the one- 
loop case under consideration is weighted by the thermal distribution either for the p meson 
or for the pion (see Hl3|). The first case is of 0(exp(— m p /T)), i.e. negligible, the second one 
therefore dominates and is expected to be of 0(T 2 /m 2 p ), due to the presence of the T = p— 
propagator in the loop (Fig. 2). 

In the following we estimate fl2.14|) in the limit m p 3> T, having in mind the mass region 
M < m p . To obtain the result analytically we make the further simplifying assumption of 
dileptons at rest (cf. (|2.5| )). Then it is straightforward to find 



i5Ujr tSE (M,q = 0) = 
e 2 J < E ') < M ~ E ') 5 '( M /2 - Bf) {E'/M - 1/4) ReIT(M - E', -p>) , (2.15) 

where \p'\ = E' . 
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When evaluating ~ReU(p ,p) on-shell, we get 



Ren(M/2,p) ~ 0{M 2 (T/m f 



(2.16) 



for \p\ = M/2, and therefore it will be neglected. The leading term of ( |2.15| ) is found after 
partial integration with respect to the variable E'\ it is proportional to the derivative 



d 
dE' 



Ren 



9; 



An 2 



T 2 
M-, 

mi 



(2.17) 



when evaluated on-shell because of the 5— function constraint in (|2.15|) , and using Boltzmann 
distributions. As expected we obtain 



i6I qfrtSE 



(M,q = 0) 



e 2 9l „ 2 T 2 
p M — 

mt 



8ir 2 4vr 



exp(-M/T) , 



(2.18) 



where we already include a factor 2 for the diagram Fig. 2 after interchanging p with p'. 

In an analogous treatment we evaluate the virtual T— dependent contributions from the 
vertex type diagrams, namely from the first two of Fig. 3. After a lengthy calculation the result 
is 



^jyBorn+ virtual 



dM 2 d A 



x 



dN 



Born 



dM 2 d 4 x 



l + f^(— ) 2 (- 
47r K m n 



19 M 2 „ 
I P 

37r 6irT 2 



n(k) 



kdk 



(k 2 - M 2 /4)' 



dN 



Born 



dM 2 d A 



x 



7 9t 



T 



1-- ^ (— ) 2 

7r Air m n 



(2.19) 



valid for m p 3> T, and for M > 0. 

In Fig. 5 we see that e.g. for T = 150 MeV the T— dependent virtual corrections ( |2.19|) 
are negative and rather large. It suggests to perform resummations, which, however, we do not 
attempt. 



3 Out of equilibrium effects 

In the early stage of heavy ion collisions, when lepton pairs are expected to be predominantly 
produced |], 0, [3|, f|, it is more likely that the pion gas is not yet in thermal and chemical 
equilibrium p5[ ^6 . Before we treat the realistic situation of an expanding gas in the next 
section, we first have to compute the production rate using non-equilibrium distributions, thus 
generalizing the calculation of sect. 2. We proceed with a tractable, but still realistic ansatz 
25] [2"6] by introducing the fugacity parameter |27J A and by using the following distributions^] 



fi(|A;o|) ~ Ae~' fc °'/ T in Boltzmann approximation (cf. (|A.6|) and ( |A.7p ). We have in mind that 
5X = A — 1 is small 5X 1 , i.e. a situation not far off equilibrium, in order to be consistent 
with the simplifying approximations summarized in Appendix A. 



2 We do not distinguish the fugacity parameter for tt^ 1 and p respectively 
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In order to estimate non-equilibrium effects even under these approximations it is not jus- 
tified to simply multiply the real emission rate ( |2.12| ) by a factor (1 + 25X) as it is in the case 
of the Born rate (cf. ( |2.7| )). Due to the structure of the (one- loop) selfenergy correction to the 
pion propagator ( |B.3| ) a more careful derivation is required. For the real emission contribution 
it first amounts to the replacement of (|2.9|) by 



i5If l2 off (q°,q) 



-xe) 



(2*0 



V 



+ ( p2) 2 + (m)2 (n(po)in 21 (p) - (1 + n(p ))iU 12 (p))\iD 21 (p'), (3.1) 

where again a corresponding term with (p <-> p') has to be added. In the second term in ( |3.1| ) 
the pinch singularity, which is not cancelled in case of non-equilibrium distributions |[28|| , is 
regularized f29| by the damping rate 7 of the pion, which we estimate in Appendix C. The 
expression for illi 2 (p) is the same as in ( |2.10|) , except that it now has to be evaluated using the 
non-equilibrium distributions ( |A.6| ), as it is the case for n(p ) in (|3.1|). 

We here discuss further details of the estimate for the dominant contribution, which is due 
to up — > 7T7* in the limit 7 — > 0, in order to exhibit our treatment, i.e. the regularization of 
the pinch singularities in the vincinity of t/u ~ (cf. Fig. 8). 

It is crucial to consider the following product of distributions, which is present in ( |3.1|) , 
(1 + n(p' )) [n(p )n(k )(l + n(p + fc„)) - (1 + n(po))(l + n(k ))n(p + fc )] , (3.2) 



adjusted for the kinematics of the process vr(p in ) + p(k p ) —>■ 7i{p ut) + 7(9), i-e. p' — > —pin, k — > 
—k p ,p — > —pout + k p (cf. Fig. 2). In relation to the definitions ( p.4|) we have to continue: 
k — > —k p ,p + — > —p out ,p- — > Pi n , which corresponds to the physical region (ii) with u > m 2 p 
in Fig. 8. The limit 7 — > is sensitive to the behaviour near t = (k p — p ou t) 2 — 0. In the 
simplifying kinematical approximation of the p - meson at rest, E p ~ m p , used in Appendix D, 

the region t ~ corresponds to positive energy po ~ t -rrz- E - > 0. 



2m D 



Thus the product of distributions 



reads in the Boltzmann approximation, 



,9/. 1 \ E P + E out 

Y(\ - l)e t 



u+m p 



(3.3) 



valid for po > and in leading order of 5X . The dominant contribution arising from 
the dilepton spectrum becomes 



to 



dN pinch 

dM 2 d A x 



7^0 



8X 



« 2 (^/47T) 7TT3 

48tt 4 M 2 \ 2m 3 p 



r 2m 2 p +M 2 



u + m p 

du e~ 2m P T 



dt+ / due'^ 7 

J2m 2 +M 2 



dt 



m 2 p M 2 



t 2 + 



9 2 P T 2 
Aire 



2 ' 



(3.4) 



where we take into account by a factor of 2 both regions denoted by (ii) in Fig. 8. The 
boundaries of the integrations are defined as in (|2.13|) , except for t <-> u. When using in (|3.4j) 
the momentum independent value for the damping 7 ( |(J . 5| ) the integrations can be explicitly 
performed. The leading term with 7 — > is proportional to I/7, indicating the pinch singularity, 



dN pinch 

dM 2 d 4 x 



- SX 



« 2 (^/4vr) 



TFT 3 



m" 



2m 2 +M 2 



24tt 3 \ 2m 3 p g 2 p /(4ire)T 



(3.5) 
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which is actually independent of the coupling g 2 . From Fig. 8 one can read off that the pinch 
near t ~ is present starting at u > m 2 + M 2 . This and (|3.3| ) explains the argument of the 
exponential function in ( |3.5| ). 
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Figure 6: Estimates of out of equilibrium effects to the dilepton rate up to 0(g 2 ) at fixed T = 150 



MeV: solid curve corresponds to /i 
to \i = 100 MeV, respectively 



; dashed to /i = 25, dashed-dotted to n = 50, and dotted 



We remark that (|3.5| ) differs numerically by only a few percent from the more exact ex- 
pression starting from ( |3.1| ) and including all the contributions due to the processes (i) - (iii) 
defined in Appendix D. 

In Fig. 6 we plot the total rate dN/dMd 4 x, including the Born term with 0(g 2 ) correction, 
for non-equilibrium cases characterized by different values of the chemical potential /i, including 
even as large values as fi = 100 MeV, which corresponds to a large value of A ~ 1.95 at T = 150 
MeV. Although here 8 A may be considered too large to justify keeping only 0(8X) terms 0, as 
done in the calculation of the curves plotted in Fig. 6, we observe that instead of an expected 
increase by a factor of A 2 ~ 4, only an effective increase of the rate by a factor 2 results, because 
of the negative contribution estimated in ( |3.5| ), when we change /i from /x = to fi = 100 MeV. 
From this fact we stress the importance of taking into account the non-trivial term (p.l|), which 



3 A different estimate in case of large A, namely replacing 7 
actually gives the same result as in Fig. 6 



A7, and Born{\ + 0(g 2 p )) 



X z Bo 



(l-O(Xg'i)) 



etc. 
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is traced back to the selfenergy diagram (Fig. 2). For completeness we state that no such extra 
terms arise from the vertex type diagrams of Fig. 3. 



4 Soft dileptons in the expanding hadronic medium 

In the previous sections we have derived the dilepton production rate originating from ir — p 
interactions up to 0{g 2 p ) in a thermal medium as a function of temperature and for fixed 
chemical potential. In the description of dilepton spectra in heavy ion collisions we have in 
addition to take into account that dilepton pairs are emitted from a rapidly expanding and 
thereby cooling medium. Hot hadronic matter produced in heavy ion collisions undergoes 
an expansion, which leads to a space-time dependence of the thermal parameters. In the 
expanding system we therefore have to integrate the rates derived in the previous sections over 
the space-time history of the collision. As a model for the expansion dynamics we assume 
the Bjorken model for 1+1 dimensional longitudinal hydrodynamical expansion f30|. In this 
model all thermodynamical parameters are only a function of the proper time, r, and do not 
depend on the rapidity variable, y. We furthermore assume, that at some initial time r ~ 1 
fm matter is formed as a thermalized hadronic gas with initial temperature To. The system 
subsequently expands and cools until it reaches the freeze-out temperature of Tf ~ 130 MeV. 
We have modeled the hadronic gas equation of state with a resonance gas where all hadrons and 
resonance states with masses up to 2.5 GeV are included. The value of the initial temperature 
T Q ~ 210 MeV has been fixed by the requirement to reproduce the multiplicity of 150 charged 
pions in the final state at central rapidity. These numbers are for S-Au collisions |J. 

With the above defined model for the expansion dynamics the space-time integration of the 
rates, first calculated in equilibrium, namely the Born rate including the virtual correction of 
0{g 2 p) ( p. 19 ), and the real processes fl2.12| ), can be performed leading to the spectrum which 



may than be compared with the one experimentally measured in heavy ion collisions. In Fig. 7 
we show the overall thermal dilepton rate from ir — p interactions including acceptance and 
kinematical cuts of the CERES experiment || . One can see in Fig. 7 that with the thermal 
source for dielectron pairs due to 7r — p interactions discussed in sect. 2 one can partially 
account for the access reported by the CERES Collaboration as measured in S-Au (and Pb- 
Au) collisions. This is particularly the case in the vicinity of the p peak where the agreement 
with the experimental data is quite satisfactory. However, the enhancement below the p peak 
and the structure of the distribution can not be explained by equilibrium production alone. 
In the kinematical window 0.25 < M < 0.45 GeV there are still almost 70% deviations of 
our theoretical curve from the experimental results. Assuming non-chemical equilibrium in the 
mesonic medium (as discussed in sect. 3) should certainly increase the dilepton rate. This is 
mostly because, in this case, there are more soft pions and rho mesons, which should produce 
more abundant soft dileptons. In Fig. 7 we show the thermal rate (long-dashed curve) assuming 
deviations from chemical equilibrium with the value of the chemical potential /i = 100 MeV. 
Indeed, we observe an increase of the rate below the rho peak. This increase, however, is rather 
modest particularly in the soft part of the spectrum. The main reason is due to the "pinch 
singular" term ( |3.5| ), which being negative reduces the contributions of the two- loop processes, 
when considered out of equilibrium. In addition the out of equilibrium distributions of pions 
and rho mesons lead to a lower initial temperature and a shorter lifetime of the thermal system. 
All these effects are the reason of small modifications of the dilepton yield. 
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Here, also additional assumptions e.g. on in-medium effects of the p meson mass and width 
eventually help to understand more completely the dilepton access measured by the CERES 
Collaboration % g 0. 



(dN e e /dr7dM)/<dN ch /d77>( lOOMeV" 1 ) 



T =1fm T f =130MeV dN ch /dy=150 



M [GeVj 
_J i L_ 



CERES 



o . o 



. 2 



. 4 



. 6 



. 8 



1 . 



1 . 2 



1 . 4 



Figure 7: Invariant-mass dielectron spectra measured by the CERES Collaboration for the 
S-Au collisions fi^f in comparison with the thermal yield from tt — p interactions. Dashed-dotted 
line describes the sum of all contributing reactions from tt^ — p° interactions. The dashed line, 
calculated by CERES, represents the expected dielectron yield from all known hadronic sources. 
The full-line is the sum of the contributions described by the dashed-dotted and dashed curves. 
The long-dashed curve includes the out of equilibrium effects 



5 Summary 

We have calculated the thermal production rate for soft dielectrons produced in a pion gas 
including relevant reactions from tt — p interactions. The calculations include the 0(g 2 p) cor- 
rections arising from the two-loop contributions to the virtual photon selfenergy. The dilepton 
rate is first calculated using equilibrium distributions. Next off equilibrium contributions to the 
rates have been estimated. Finally, dilepton production in an expanding hadronic medium has 
been calculated applying the Bjorken model. No medium effects on he p resonance parameters 
are included in the analysis. The access reported by the CERES Collaboration in the low mass 
range below the p peak is diminished, although not yet completely described by the model. 
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Appendix A: Thermal propagators 

In order to perform the perturbative calculations presented in this paper we use the closed- 
time-path formalism of thermal field theory |18, [l^, ^(J. The real-time 2x2 matrix (scalar) 
propagators, 

Du(fc) = (l + n(k ))A R {k)-n{k )A A {k), 
D 12 (k) = n(k )(A R (k)-A A (k)), 
D 21 (k) = (l + n(k ))(A R (k)-A A (k)), 

D 22 (k) = n(k )A R (k)-(l+n(k ))A A (k), (A.l) 
are expressed in terms of the retarded and advanced propagators, 

where k 2 = k 2 , — k 2 and m is the T = mass. For equilibrium conditions the Bose-Einstein 
distribution function n(ko) depends on temperature T, 

and satisfies 

n(k )+n(-k ) + 1 = 0. (A.4) 
For our numerical estimates we are mainly using the Boltzmann approximation 

n(|fc |)~exp(-N/T). (A.5) 



Estimating non-equilibrium effects we follow the approximations described in detail in ||19|| . 



and more recently in pi| . It amounts to replace the thermal n(ko), ( |A.3| ), by their non- 
equilibrium counter parts, i.e. in general by Wigner distributions n(k ,X). This corresponds 
in terms of the cumulant expansion to approximate non-equilibrium correlations by the second 
cumulant only. As a further approximation we suppress the possible dependence on the center 
of mass coordinate X, essentially assuming a homogenous and isotropic medium. 

Practically for our estimates we take (cf. (|A.4j )), 

n(A; ,X)^ ko> °> (A.6) 

\ -(l+n(N) k <0, y ' 

with 

= ePo|-, 1 )/T_ 1 = e | fc ol/T_ A ' ( A ' 7 ) 

by introducing the fugacity parameter p7| A = e M / T , which is assumed to be energy indepen- 
dent. Obviously A ^ 1 in case of (chemical) non-equilibrium. Consequently deviations from 
equilibrium are approximated in terms of the equilibrium distribution ( |A.3| ) by 

071 

—— = n(k , X)/n(k ) -1~8\ e(k ) (1 + n(k )) (A.8) 
n{Ko) 
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in leading order of SX = A — 1. 

In the expressions for A^( A )(/c) the usual limit e — > is considered, whenever this limit is 
defined, e.g. 

A R (k) - A A (k) = -2me(k )5(k 2 - m 2 ), (A.9) 

where e(k ) = 9(k ) — 9(—k ) is the sign function. 

However, out of equilibrium pinch singularities appear in this limit e — > PB| , e.g. in 
products like An(k)AA_(k) (see Appendix B). We treat this situation by following the conjecture 
by Altherr [JIJ keeping e non-vanishing: it is identified by the damping width 7 > 
(Appendix C). 



Appendix B: One-loop selfenergy correction to the pion propagator 

Here we comment (cf. Fig. 2) on the selfenergy correction of the pion propagator and the pos- 



sible appearence and presence of pinch singularities |18|, |28| . Denoting the selfenergy insertion 
by the matrix n a f,(p) we obtain the improved (12)-propagator using the Dyson equation, 

iD 12 {p) -> iD 12 (p) + J2 iD la (p)(-iIl ab (p))iD b2 (p) 

a,b=l,2 

= iD 12 (p)+i5D 12 (p). 
With the propagators specified in the Appendix A, and with the relations 



u n (p) = -n* 2 ( P ), imn n = -(n 12 + n 21 ), 

also valid out of equilibrium, we obtain 

5D 12 (p) = n(p )(A 2 ,(p)-A 2 4 (p))Ren 11 (p) 

+ ln(p )(A 2 R (p) + A 2 A (p))[Tl 12 (p)-Tl 21 



(B.l) 
(B.2) 



+ A fl (p)A A (p)[n(po)n 21 (p) - (1 + n(p ))U 12 (p)]. (B.3) 

Here the ill-defined product Ar(p)Aa(p) appears, giving rise to possible unpleasant pinch 
singularities JES), In case of equilibrium, however, it is well known that the last term in 
(B.3) vanishes due to the detailed balance condition, 

n 2 i(p) = e Po/T Tl 12 (p). (B.4) 

As a consequence (B.3) simplifies to 



5DH(p) = 2irie(p )n(po) 



S'(p 2 )ReU(p) + -P — Imn(p) 



where 



ReU(p) = Re nn (p), Imn 



7T \P , 

1 JUuip) 

2 n{po) 



(B.5) 



(B.6) 



expressed in terms of the equilibrium distribution n(po), and in terms of the (retarded) self- 
energy H(po + iepo,p) f2~3fl . 5' denotes the derivative of the 5-function, and P the principal 
value. 

However, for nonequilibrium distributions ( |A.6| ), for which obviously (l+n)/n 7^ exp(p /T), 
the last term in (B.3) does not cancel! In order to regularize the product Ar(p)Aa(p) |29j we 
take into account the non-vanishing (on-shell) damping rate of the pion, which we estimate in 
the following Appendix C. 



16 



Appendix C: Estimate of the damping rate of the pion 

The pion damping rate 7 is determined by the "pole" (in the lower energy half-plane) of the 
retarded propagator. In the following we not only neglect the pion mass, but also thermal 
corrections due to the real part of the (retarded) pion selfenergy, such that 

lmU(p ,p) 

7— > l^- 1 ) 

Po 

evaluated on-shell p 2 = for positive pion energy p . In the one- loop approximation the 
absorptive part of the pion selfenergy is given by 

9p f d^k , . / k a k T 



lmU{P0 ^ = 2^[)I ^^P + kr[-9,r + ^j(2p + kyD 12 ( P + k)D 21 (k), (C.2) 

where the momentum labels of Fig. 2 are used. 

The temperature dependent non- vanishing part for p 2 = p\ — p 2 = 0, becomes (cf. the 
derivation in fll3[| ), 



r d 3 k 1 

lmU(p =p,\p\=p) = -g 2 p m 2 p n J — — [n(K) ~ n{E p )\ 5(p - E p + E«), (C.3) 



(2tt)3 2^2^ 



with the corresponding energies are E p = J k 2 + m 2 p and E n = \p + k\. 

The dominant contribution comes from the pion's thermal distribution n(E n ). In the Boltz- 
mann approximation it leads to 

9 9 9 9m 2 

(7 777, POO Q Tfl _[ m p 

lmU( Po = p,p)^-^ j^n(E„)dE^-^^e-^. (C.4) 

Imll vanishes for p — 0, and it has its maximum near p ~ We note that ( |C.4| ) gives a 
positive damping rate 7, as required. In order not to overestimate the contributions arising 
from the A^A^ terms we take for the numerical estimates a momentum independent value, 
namely the one at the maximum, 

9 2 1 

p 7 ~ — Imll (po ~ p,p) ~ — ^--T 2 . (C.5) 

An e 

Appendix D: Kinematic boundaries for the kinetic processes 

The dilepton thermal rate receives contributions from the two-body processes tttt — > p7* and 
up — > 7T7*, and from p — > 7T7T7*, for which we evaluate in the following the phase-space integral 
( 2. 11|) including its kinematic constraints. The calculation is simplified using (i) Boltzmann 



distributions and (ii) the fact that m p /T ^> 1: the p-meson is kept at rest in the medium. This 
allows to integrate the p's phase-space by 



and to replace elsewhere E p = J k 2 + m 2 ~ m p . 
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The remaining integrations are carried out as we explicitly show in the following for the 
channel ir(p + ) + n(p_) — > p(k) + r y*(q). We only consider the case of dilepton masses M < m p . 

The phase-space integral defined by 
is approximated by 



with q ~ y\p+ + p~\ 2 + M 2 . The angular integrations in I can be done. 




Figure 8: Mandelstam plot: physical regions for the processes (i) nil — > fry* , (ii) Tip — > 717* 
and (Hi) p — > 7T7T7* 

Next we introduce the Mandelstam variables by 

s = Go + +p_) 2 = (/c + g) 2 , 
t=(p + -A;) 2 = (g-p_) 2 , 

M =(p__A;) 2 = (q-p + )\ (D.4) 
with s + t + u = m 2 + M 2 . Approximating the pion energies by 



E + ~ (m 2 - *)/2m p , £_ ~ (m 2 - «)/2m p , (D.5) 



"p V/ -"»pj — - — V""p 

gives the final result for the phase-space integral, 



I~^—I p I dtdu 6{tu — m 2 M 2 ) 6(s — (m p + M) 2 ). (D.6) 
8mi J p 
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The other channels only differ in their kinematical boundaries, 

(ii) up -> tt 7 * : 6(-s) 6{m 2 p M 2 - tu), 

(m)p->7T7r7* : 0(s) 6{tu - m 2 p M 2 ), (D.7) 

which are also summarized in Fig. 8. The processes imp — > 7* and ir — > 7rp7* do not contribute 
due to the nonvanishing masses m p and M. 

We finally remark that the poles at t — and/or tt = are present within the kinematical 
domain for the channel up — > 717*, whereas for 7r7r — > p7* they are present just at the kinematical 
boundary when s — > 00. 
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